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The microscopic mechanism of thermal dissipation in quantum turbulence is numerically studied by 
solving the coupled system involving the Gross-Pitaevskii equation and the Bogoliubov-de Gennes 
equation. At low temperatures, the obtained dissipation does not work at scales greater than the 
vortex core size. However, as the temperature increases, dissipation works at large scales and it 
affects the vortex dynamics. We successfully obtain the mutual friction coefficients of the vortex in 
dilute Bose-Einstein condensates dynamics as functions of temperature. 
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The physics of quantum turbulence (QT), comprising 
tangled quantized vortices, is one of the most important 
research topics in low-temperature physics Stimu- 
lated by recent experiments on both superfluid 4 He and 
superfluid 3 He where a few similarities were observed be- 
tween quantum and classical turbulences El ED > stud- 
ies on QT have entered a new stage where one of the 
main motivations is to investigate the relationship be- 
tween quantum and classical turbulences. 

With this motivation, we theoretically investigated QT 
in our previous studies by numerically solving the Gross- 
Pitaevskii (GP) equation |SJ|S]. Since the GP equation is 
applicable to compressible quantum fluids, compressible 
excitations of wavelengths smaller than the vortex core 
size, which affect the vortex dynamics in QT, are emitted 
during vortex reconnections, the disappearance of small 
vortices, or by high frequency Kelvin waves. As a result, 
the intrinsic behavior of quantized vortices in QT, such 
as the Richardson cascade, is hindered by the compress- 
ible short-wavelength excitations. To eliminate these ex- 
citations, we introduced a phenomenological dissipation 
term that works only at scales smaller than the vortex 
core size and successfully obtained the Kolmogorov law in 
QT, which is one of the most important statistical laws in 
classical turbulence [7|. However, the microscopic origin 
of dissipation and the realistic nature of the introduced 
dissipation term in quantum fluids were still unknown. 

In the research field of superfluid 4 He, the dissipation 
mechanism in quantized vortices has been studied for ap- 
proximately 50 years with the consequent development 
of the concept of mutual friction. The hydrodynamics 
of superfluid 4 He is usually described using the two-fluid 
model; in this model, the system consists of an inviscid 
superfluid and viscous normal fluid, and dissipation in 
quantized vortices is caused by the mutual friction be- 
tween the vortices and normal fluid jS] . Although mutual 
friction is very important in QT EH HH E2 EH , no 
microscopic theory exists that can reproduce the mutual 
friction coefficients as functions of temperature. This is 
because the superfluid 4 He system is strongly correlated 
and extremely difficult to study microscopically by using, 



for example, the quantum field theory. 

Quantized vortices were also discovered in systems of 
atomic Bose-Einstein condensates (BECs) ^1 ■ Since 
the systems are different from superfluid 4 He, they can 
be easily studied by using the GP equation for the con- 
densate and the Bogoliubov-de Gennes (BdG) equation 
for the excitations from the condensate; these equations 
should also prove useful in revealing the microscopic dis- 
sipation mechanism in quantized vortices. However, in 
the research field of atomic BECs, no concept exists that 
explains dissipation or mutual friction of quantized vor- 
tices. Some groups E3E1E3 nave discussed dissipation 
in atomic BECs; however, they have not directly consid- 
ered the dissipation mechanism in quantized vortices. In 
this work, we are the first to study the dissipation mech- 
anism in quantized vortices in dilute BECs. With this 
study, we cannot only justify the phenomenological dis- 
sipation term introduced in our previous work but also 
define the mutual friction in a dilute BEC by considering 
its relationship with the mutual friction in superfluid 4 Hc. 

It would be fairly reasonable to assume that dissipa- 
tion in quantum fluids, including the quantized vortices 
is caused by the interaction between the condensate and 
its excitations. The dynamics of the excitations can be 
described by the BdG equation; therefore, we numerically 
solve the time development of the coupled system involv- 
ing the GP and BdG equations for a quantum fluid with 
some quantized vortices. Our results reveal that dissi- 
pation obtained by the coupled equations can work only 
at small scales and reproduces the phenomenological dis- 
sipation term introduced at very low temperatures, as 
mentioned in our previous studies. As the temperature 
increases, however, dissipation works at larger scales and 
it affects the vortex dynamics; this is qualitatively similar 
to the mutual friction of quantized vortices in superfluid 
4 He. Moreover, our model is limited to a dilute Bose 
gas. Thus, we further investigate dissipation as a model 
for mutual friction in dilute BECs by calculating the dy- 
namics of one straight vortex, and we successfully obtain 
the friction coefficients as functions of temperature. It is 
impossible to directly compare the obtained friction coef- 
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ficients with those of superfluid 4 He because both the GP 
and BdG equations are applicable only to dilute BECs. 
However, this study can be regarded as the first study 
elucidating the mutual friction of quantized vortices in 
dilute BECs. 

To consider a quantum fluid as a Bose-Einstein con- 
densed system, we start with the many-body Hamilto- 
nian, 



H 



dx * 



Its dynamics can be described by 
.9* 



Of 



[-V 2 - fi + 0*t#]tf. 



(1) 



(2) 



Here, if?(x,t) is the boson field operator, /x the chemi- 
cal potential, and g the coupling constant. In the Bose- 
Einstein condensed system, the field operator if?(x,t) can 
be represented in terms of the mean-field ansatz |20j , 



* = $ + x + C, 



(3) 



which conserves the particle number. Here, we define 
the terms appearing in Eq. ((3J): the macroscopic wave- 
function is represented as Q(x, t) = O(y/No/V), the first- 
ordered excitations as x{x,t) — 0(1/ W), the higher- 
ordered excitations as ((x,t) = O(l/y/N V), the num- 
ber of condensate particles as Nq, and the volume of the 
system as V. Substituting Eq. © into Eq. (0) and ne- 
glecting the higher-ordered excitations ((x,t), we obtain 
the GP equation as 

= hV 2 — H + 5 (|$| 2 + 2<x t x))]4> + g(x 2 W (4) 
and the BdG equation |^ [22] as 



/i + 2 5 |$| 2 ]x + 5* 2 X t 



(5) 



When the macroscopic wave-function is expressed as 
$(x,t) = f(x,t)e i ^ x ' t \ f(x,t) 2 is considered to be the 
condensate density, and the superfluid velocity v(x,t) is 
given by v(x,i) = 2V(j>(x,t). The vorticity u:(x,t) — 
V x v(x,t) vanishes everywhere in a single connected re- 
gion of the fluid and thus any rotational flow is carried 
only by quantized vortices. In the core of each vortex, 
&(x, t) vanishes; therefore, the circulation § v ■ ds around 
the core is quantized by An. The vortex core size is given 
by the healing length £ = l/fy/g. 

The GP equation (@J| can be expressed as 
id$(x,t)/dt = H GP $(x,t). The Hamiltonian H GP 
of the GP equation has the following imaginary term: 



Im[iJGp] = — 7 = fflnr 



(6) 



This defines the dissipation 7(2;, t) of the condensate 
caused by the interaction with the noncondensed par- 
ticles. We can quantitatively calculate the dissipation 



7(2:, i) from Eq. © by numerically solving the coupled 
system involving the GP equation and BdG equation 

©■ 

To solve the BdG equation (J5J, we use the Bogoliubov 
transformation; therefore, 
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(7) 



where uj(x, t) and vj(x, t) are the Bogoliubov coefficients, 
and &j and at are the annihilation and creation opera- 
tors of a quasiparticle, respectively. We assume that the 
quasiparticles are coupled with a heat bath at tempera- 
ture T; further, we use the local equilibrium assumption 



{a)a 3 ) = Nj = 



1 



exp{Ej/T) - 1 



(8) 



with the excitation spectrum E 3 of quasiparticles. The 
local equilibrium assumption © yields the energy flow in 
the system, as shown in Fig. ^ When excitations such 
as vortices or sound waves are formed in the condensate, 
their energy is transferred to quasiparticles and finally 
dissipated to the heat bath. By using Eqs. J7J) and JSJ, we 
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Figure 1 : Image of the energy flow in a system of quan- 
tum fluid coupled with a heat bath. 

can deduce the final form of the coupled system involving 
the GP and BdG equations: 



d<f> 
l dt 

.dvj 



E; = 



[-V 2 - v + 5 (|$| 2 + 2n )]$ + «?m e $*, (9a) 

[-V 2 - M + 2 5 |$| 2 ] % - g®\ = A 3 , (9b) 

- [-V 2 -n+ 2g\$\ 2 ]vj + g$* 2 u 3 = B h (9c) 

^[KfiV,- + Kf(^ + 1)], (9d) 
3 

-Y^UjV^Ni + l)], (9e) 
3 

J dx Re • ,•;//,'. (9f) 



We numerically solve the coupled equations (|5J. We 
use a pseudospectral method [23] in space with periodic 
boundary conditions in a box with a spatial resolution 
of 32 3 grid points. With regard to the numerical pa- 
rameters, assuming g = 1, we use a spatial resolution of 
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Am = 0.125 and V — 4 3 , where the length scale is nor- 
malized by £. The numerical time evolution is given by 
the Runge-Kutta-Gill method [23] with a temporal reso- 
lution of At — 1 x 1CP 4 . We begin with a macroscopic 
wave-function $(x,t — 0) that includes several randomly 
placed vortices, as shown in Fig. [21 (a), and with the 
uniform excitations x(x, t — 0) given by, 



Uj(x,t = 0) = e 1 ^ 



2V Ej 



Vj(x, t — 0) = e 



21/ E. 



(10) 



where fc^ = 2irj /-^V, where j is an integer. Subsequently, 
we calculate the dissipation term j(x,t) at t = 1. Figure 
121 (b) shows the dependence of the Fourier-transformed 
dissipation y(k,t) on the wavenumber k at several tem- 
peratures: T = 0.01T C , 0.1T C , and 0.5T C , where T c = 
4-7r/{£(3/2)} 2 / 3 is the critical temperature for the Bose 
condensation of free bosons. At a low temperature of 
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Figure 2: (a) An example of the configurations of quan- 
tized vortices at t = 0. (b) Wave number dependence 
of the Fourier-transformed dissipation term j(k,t = 1), 
which is obtained by performing an ensemble average of 
25 initial states. 

T = 0.01T C , dissipation works only at wave numbers 
greater than 27r/£, which is consistent with the dissipa- 
tion term j(k) = r yo9(k — 27r/£) with a step function 9 
introduced in our previous studies |SJ E] . From this re- 
sult, we find that only short-wavelength excitations emit- 
ted during the vortex reconnections, the disappearance of 
small vortices, or by high frequency Kelvin waves get dis- 
sipated at scales smaller than the vortex core size. This 
phenomenon has also been confirmed by another simu- 
lation that starts from a state without vortices and ob- 
tains smaller and less fc-dependent values of "f(k,t). On 
the other hand, as the temperature increases, dissipation 
works at small wave numbers as well. Since dissipation at 
small waven umbers dissipates vortices at scales greater 
than the vortex core size £, the vortex dynamics are di- 
rectly affected by this dissipation. Hence, we can expect 
an effect of "f(x, t) similar to that of mutual friction in 
superfluid 4 He. 



Figures |3| (a) and[3J(b) show the vortex configurations 
at t = 1 and for T = 0.01T C and T = 0.1T C respec- 
tively, starting from the vortex configuration shown in 
Fig. |21 (a) . Because dissipation at large scales affects both 
small vortices and short-wavelength Kelvin waves along 
the vortices, the vortex configurations clearly show the 
difference in the effects on vortex dynamics at different 
temperatures. We can observe fewer vortices and fewer 
Kelvin waves in Fig. 01(b) than in Fig. 01(a), which obvi- 
ously conforms to the effect of dissipation at large scales. 
By using the vortex-filament model, Tsubota et al. in- 
vestigated the difference in the vortex configurations in 
QT at different temperatures and observed fewer Kelvin 
waves at higher temperatures |24j . Their observation is 
qualitatively consistent with our result. 
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Figure 3: Configurations of quantized vortices at t = 1 
at T = 0.01T C (a) and T = 0.1T C (b), starting from the 
configuration shown in Fig. [21(a). 

Next, we attempt to calculate the coefficients of mu- 
tual friction as functions of temperature. When one 
straight vortex along the z-axis is placed under the ve- 
locity field v c = [v c , 0, 0], the dynamics of the vortex po- 
sition s(t) = [s x (t), s y (t),0] can be described by mutual 
friction as s(t) — [a'v c ,av c ,0] miHl- Here, a and a' are 
the coefficients that characterize the amplitude of mutual 
friction. The solution becomes 



s(t) = [s x (0)+a'v e t,s y (p)+av e t,Q] 



(11) 



Starting from the state with one straight quantized vor- 
tex, we numerically solve the coupled equation involving 
the GP and BdG equations under the velocity field 



9$ 



[-V 2 — fi + g(\<P\ 2 + 2n c ) + iv c ■ V]$ 
[-V 2 -/x 



(12a) 

2#|$| 2 + iv c ■ V]uj ~g^ 2 Vj, (12b) 
2g\<&\ 2 + iv c ■ Vjvj 



(12c) 



for the case of v c = 0.1. We can calculate a and a' by 
comparing the position of the vortex in the numerical 
simulation using Eq. (|llfl . Figure^ shows the tempera- 
ture dependence of a and a' . The effect of the dissipation 



3 



term 7(2;, t) on the vortex dynamics at large scales and its 
monotonic increase with temperature are evident, which 
is qualitatively consistent with mutual friction in super- 
fluid 4 He at temperatures much lower than the superfluid 
critical temperature. This result develops the first esti- 
mation of mutual friction of quantized vortices in dilute 
BECs. The temperature dependence of a and a' needs to 
be experimentally observed, and it may become a stan- 
dard scale for measuring the temperature in atomic BECs 
with quantized vortices. 
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Figure 4: Temperature dependence of the coefficients a 
and a'. Plots represent the numerical results and lines 
indicate some fitting. 

This work allows quasiparticles to move via Eq. J5J, 
but we do not discuss the motion here. With regard to 
turbulence in superfluid 4 Hc, Vinen predicted that at fi- 
nite temperatures, superfluid and normal fluid are likely 
to be coupled together at large scales due to mutual fric- 
tion and thus behave similarly to the turbulence in a one- 
component fluid |25j . In our simulation, we can also ex- 
pect a similar coupled turbulence in which the dynamics 
of the quasiparticles is strongly coupled with that of the 
condensate with both the dynamics becoming comparable 
at large scales. This physics will be reported shortly. 

In conclusion, we investigated the dissipation mech- 
anism in quantized vortices by numerically solving the 
coupled equation involving the GP and BdG equations. 
At low temperatures, dissipation works at scales smaller 
than the vortex core size, which is consistent with the 
phenomenological dissipation term introduced in our pre- 
vious studies. As the temperature increases, dissipation 
can work at large scales as well, and it directly affects the 
vortex dynamics. We successfully related this effect to 
the mutual friction in superfluid 4 He by calculating the 
mutual friction coefficients as functions of temperature. 
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